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In order to gain a deep insight into /-electron properties of filled skutterudite compounds from a micro- 
scopic viewpoint, we investigate the multiorbital Anderson model including Coulomb interactions, spin- 
orbit coupling, and crystalline electric field effect. First we examine the local /-electron state in detail in 
comparison with the results of LS and j-j coupling schemes. For each case of n=l~13, where n is the 
number of / electrons per rare-earth ion, the model is analyzed by using the numerical renormalization 
group (NRG) method to evaluate magnetic susceptibility and entropy of / electron. In particular, for the 
f 2 -electron system corresponding to the Pr-based filled skutterudite, it is found that magnetic fluctuations 
significantly remain at low temperatures, even when the ground state is "i singlet, if r 4 triplet is the 
excited state with small excitation energy. In order to make further step to construct a simplified model 
which can be treated even in a periodic system, we also analyze the Anderson model constructed based on 
the j-j coupling scheme by using the NRG method. It is clearly observed that the magnetic properties are 
quite similar to those of the original Anderson model. Then, we construct an orbital degenerate Hubbard 
model based on the j-j coupling scheme to investigate the mechanism of superconductivity of filled skut- 
terudites. In the 2-site model, we carefully evaluate the superconducting pair susceptibility for the case of 
n=2 and find that the susceptibility for off-site Cooper pair is clearly enhanced only in a transition region 
in which the singlet and triplet ground states are interchanged. We envision a scenario that unconventional 
superconductivity induced by magnetic fluctuations may occur in the / 2 -electron system with T\ ground 
state such as Pr-based filled skutterudite compounds. 

KEYWORDS: Filled skutterudites, Multiorbital Anderson model, Numerical renormalization group 
method, Magnetic fluctuations, Orbital fluctuations, Unconventional superconductivity 



1. Introduction 

Recently /-electron compounds with filled skutteru- 
dite structure, expressed as RT4X12 (R=rare-earth ion, 
T=transition metal ion, and X=pnictogen), have attracted 
much attention in the research field of condensed matter 
physics. 1 Since the filled skutterudite compounds exhibit 
huge thermopower, they have been focused as thermoelec- 
tric materials, mainly for the purpose of industrial applica- 
tion. However, after the discovery of heavy fermion behavior 
and superconductivity in PrOs4Sbi2, 2,3 the filled skutterudite 
compound has been an important target of active investiga- 
tions also from the viewpoint of basic material science, since 
this is the first Pr-based / 2 -electron system that shows heavy- 
fermion superconductivity. 

From experimental results on specific heat, 4 magnetiza- 
tion, 5 and neutron scattering, 6 ' 7 it has been recently confirmed 

that the ground state of PrOs4Sbi2 is Ti singlet, while the 

(2) 

first excited state is T\ triplet with very small excitation en- 
ergy less than 10K. Note that the filled skutterudite has cu- 
bic symmetry described by not Oh point group. 8 In such 
a material with non-magnetic Ti singlet ground state, con- 
ventional superconductivity is naively expected. However, the 
superconductivity of PrOs4Sbi2 is suggested to be unconven- 
tional from NMR experiment, 9 since there is no coherence 
peak in 1/Ti, where T\ is nuclear spin-lattice relaxation time, 
just below the superconducting transition temperature T c . 

Here one may expect that PrOs4Sbi2 belongs to the group 
of unconventional d-wave superconductors, as has been fre- 
quently observed in strongly correlated electron systems in- 
cluding Ce- and U-based heavy fermion superconductors. 10 



However, 1/Ti of PrOs4Sbi2 does not obey the so-called T 3 
behavior, 9 characteristic of the line-node gap function in the 
d-wave pairing. Rather, the isotropic superconducting gap is 
consistent with the experimental results, in which 1/T\T ex- 
ponentially decreases even crossing T c . In addition, the su- 
perconductivity of filled skutterudite compounds exhibit sev- 
eral exotic features, for instance, point-node behavior in the 
gap function, 11, 12 multiple superconducting phases, 11 and the 
breaking of the time-reversal symmetry detected by /iSR ex- 
periment. 13 Then, it has been highly requested to elucidate 
theoretically the mechanism of such exotic superconductivity 
ofPrOs 4 Sbi 2 . 

Let us turn our attention to Pr-based filled skutterudites in- 
cluding other transition metal ions and pnictogens. It has been 
found that the ground-state properties are changed sensitively 
depending on T and X in PrT4Xi2. As mentioned above, 
PrOs4Sbi2 exhibits unconventional superconductivity, while 
PrRu4Sbi2 is confirmed to be conventional s-wave supercon- 
ductor from Sb-NQR experiment. 14 Interestingly enough, due 
to the substitution of transition metal ion, superconducting 
properties in Pr(Osi_ x Ru x )4Sbi2 are changed from uncon- 
ventional to conventional ones. 15 When pnictogens are sub- 
stituted, the ground-state nature is also drastically changed. 
For instance, PrOs4Pi2 is non-magnetic metal, 16 PrRu4Pi2 
exhibits metal-insulator transition, 16, 17 and PrFe4Pi2 shows 
exotic quadrupolar ordering. 18 It is also an important task im- 
posed on theoretical investigations to clarify what is a key is- 
sue to control the electronic properties of Pr-based filled skut- 
terudites. 

In addition to Pr-based compounds, numerous kinds of 
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filled skutterudite materials have been synthesized by ex- 
pert techniques. Those compounds actually exhibit varieties 
of electronic properties: La-based filled skutterudite mate- 
rials are known to be conventional BCS superconductors. 9 
Ce-based filled skutterudites are Kondo semi-conductors with 
the gap up to thousand Kelvins. For the skutterudites includ- 
ing rare-earth ions other than La, Ce, and Pr, the ferromag- 
netic ground state has been frequently observed, for instance, 
in RFe 4 P 12 for R=Nd, Sm, Eu, Gd, Tb, Dy, and Ho. How- 
ever, in some Sm-, Gd-, and Tb-based skutterudites such as 
SmOs 4 Sbi 2 , 19 GdRu 4 Pi2, 20 and TbRu 4 Pi 2 , 20 antiferromag- 
netic ground states have been also found. It is another char- 
acteristic issue of filled skutterudites that the magnetic prop- 
erties are changed depending on the number of /-electrons n 
per rare-earth ion. 

Even in the above brief survey of the properties of filled 
skutterudite compounds, we can easily understand that this 
group of materials exhibits richness, as revealed by exper- 
imental effort. However, theoretical investigations on mag- 
netism and superconductivity for /"-electron systems with 
n> 1 have been almost limited to the phenomenological level 
based on the Ginzburg-Landau theory, mainly due to com- 
plexity in the many-body problem originating from the com- 
petition among strong spin-orbit coupling, Coulomb interac- 
tions, and the effect of crystalline electric field (CEF). For the 
purpose to include such interactions, an LS coupling scheme 
has been widely used, but it is not possible to apply stan- 
dard quantum-field theoretical technique in the LS coupling 
scheme, since Wick's theorem does not hold. 

In order to overcome such a difficulty, it has been proposed 
to construct a microscopic model for /-electron systems from 
a j-j coupling scheme. 21 Since individual /-electron states 
are clearly defined in the j-j coupling scheme, it is conve- 
nient for the inclusion of many-body effects using standard 
theoretical techniques. In fact, based on such a microscopic 
model, the following points have been clarified: Key role of 
orbital degree of freedom for superconductivity of CeTIn5 
and PuTGa 5 with transition metal ion T, 21 " 23 an odd-parity 
triplet pair induced by Hund's rule interaction, 24 the com- 
plex magnetic structure of UTGa5 25 and NpTGa5 26 based on 
a spin-orbital model, and microscopic origin of octupole or- 
dering of Np02- 27 

When we attempt to apply the microscopic model based 
on the j-j coupling scheme to the filled skutterudite com- 
pound, unfortunately, there are several problems which should 
be clarified. First of all, the application of the j-j coupling 
scheme to the rare-earth materials seems questionable except 
for cerium compounds, since the Coulomb interactions are 
larger than the spin-orbit coupling in 4/-electron systems. 
Second, it is still unclear how to include the CEF effect for 
/"-electron systems in the j-j coupling scheme, in which the 
CEF effect is just a one-electron potential determined for the 
f 1 -electron system. 21 It is necessary to show explicitly how 
the j-j coupling scheme reproduces the electronic properties 
of /"-electron systems. 

In this paper, first we analyze in detail the /-electron term 
including spin-orbit coupling, Coulomb interactions, and CEF 
effect. In particular, we carefully compare the local /-electron 
states with those of the LS and j-j coupling schemes. We 
also show the CEF energy levels for n=l~13 under the CEF 
potential of Th point group. By further adding the conduc- 



tion electron and hybridization terms, we obtain the Anderson 
model and analyze it numerically. Then, the Curie-law behav- 
ior in magnetic susceptibility is found even in the / 2 -electron 
system, when 1/ singlet is the local ground state and 
triplet is the excited state with small excitation energy. The 
same tendency is observed in the Anderson model constructed 
from the j-j coupling scheme. Then, we further construct the 
Hubbard-like model based on the j-j coupling scheme and 
evaluate the superconducting pair susceptibility. The numeri- 
cal results on the model suggest that anisotropic Cooper-pair 

(2) 

may appear in a region in which Ti singlet and T y 4 ' triplet 
states are interchanged. 

The organization of this paper is as follows. In Sec. 2, we 
discuss in detail the properties of local /-electron term in 
comparison with those of the LS and j-j coupling schemes. 
Then, we show the energy levels by changing the /-electron 
number n between n=\ and 13 for the fixed CEF parame- 
ters of the Th point group. In Sec. 3, the Anderson model 
is introduced for filled skutterudites and it is analyzed by 
using the numerical renormalization group (NRG) method. 
The results for magnetic susceptibility and entropy are shown 
for n=l~13. In order to reduce effectively the number of /- 
orbitals, we also introduce the Anderson model in the j-j cou- 
pling scheme by focusing on the case of n=2, i.e., the case of 
Pr-based filled skutterudites. The model is analyzed by using 
the NRG method and it is shown that the essential point of the 
results in the original Anderson model can be correctly cap- 
tured even in the model based on the j-j coupling scheme. In 
Sec. 4, we construct an orbital degenerate Hubbard model ap- 
propriate for filled skutterudites. By evaluating carefully the 
superconducting pair susceptibility of the model, we clearly 
show that the pair susceptibility is enhanced only in a narrow 
region in which the ground state is changed between singlet 
and triplet states. In Sec. 5, we will summarize the paper. 

2. Local f-electron state 

In this section, let us examine in detail the properties of 
local /-electron states in comparison with those obtained in 
the LS and j-j coupling schemes. 

2.1 Local f -electron Hamiltonian 

In general, the local /-electron term is composed of three 
parts as 

Hf = H c + H so + H C ef, (1) 
where Hq is the Coulomb interaction term, written as 



,m2,m3,m4 



mi~rri4 <7\ ,<T2 



^ frn\<j\ fm2<72 f rn 3&2 fm^a i • 



(2) 



Here f ma is the annihilation operator for /-electron with 
spin a and angular momentum m(=— 3,- • • ,3) and cr=+l (— 1) 
for up (down) spin. The Coulomb integral I mi ,m 2 .m3,m4 is 
known to be written in the form of 



^2F k c k (m 1 ,m4)ck(m 2 ,m 3 ), (3) 



k=0 



where the sum on k includes only even values (k=0, 2, 4, and 
6), F k is the Slater-Condon parameter including the complex 
integral of the radial function, 29 30 and is the Gaunt co- 
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efficient, 31 - 32 which is tabulated in the standard textbooks of 
quantum mechanics. 33 It is convenient to express the Slater- 
Condon parameters as 

5° = A+ 15C + 95/7, 
F 2 = 225(5 - 6C/7 + D/42), 
F A = 1089(5(7/7+5/77), 
F 6 = (429/5) 2 • (-D/462), 

where A, 5, C, and D are the Racah parameters. 32 
The spin-orbit coupling term, H so , is given by 



(4) 



X 



t' fma /) 



mcrJra 1 a 1 ; 



(5) 



where A is the spin-orbit interaction and the matrix elements 
are expressed by 

Cm+l,l,m,T = V 12 ~ m ( m + i)/ 2 ' 



(6) 



Cm-x,T,m,4. = V 12 - rMm - l)/2, 

and zero for other cases. 

The CEF term Hqef is given by 

-HCEF = 22 B m ,m' fmafm'o, 



(7) 



where 5 m , m < is determined from the table of Hutchings for 
angular momentum £=3, 34 since we are now considering the 
potential for / electron. For filled skutterudites with Th sym- 
metry, 8 B m rn i is expressed by using three CEF parameters 

5 40 , 5 60 , and B 62 as 

53,3 = -B_ 3i _3 = 1805 40 + 180560, 

B 2 ,2 = B-2,-2 = -4205 40 - 10805 60 , 
B hl = = 605 40 + 2 7 005 60 , 

5 . = 3605 40 - 36005 60 , 

Bs,-i = 5-3.1 = 60a/15(5 40 - 21560), (8) 
5 2 ,_ 2 = 3005 40 + 75605 60 , 

5 3 .1 = 5_ 3 _i = 24a/155 62 , 
5 2 ,o = 5_ 2 ,o = -48yi55 62 , 
51,-1 = — 5 3: _3 = 3605 6 2- 



Note the relation of 5 m m /= 
notation, we define 



: 5 m ',m- Following the traditional 



5 40 = Wx/F(4), 

B m = W(l-\x\)/F(6), 

5 62 = Wy/F\6), 



(9) 



where x and y specify the CEF scheme for Th point group, 
while W determines an energy scale for the CEF potential. 
Although 5(4), 5(6), and 5* (6) have not been determined 
uniquely, 35 in Eq. (7) we choose 5(4)=15, 5(6)=180, and 
5 t (6)=24for£=3. 34 

We note that the CEF potential is originally given by the 
sum of electrostatic energy from the ligand ions at the posi- 
tion of /-electron ion, leading to the one-electron potential 
acting on the charge distribution of /-orbitals, as expressed 
by Eq. (7). Thus, in principle, it is not necessary to change 
the CEF potential depending on the /-electron number. As 
we will see in the following subsections, the CEF schemes 
for n=1^13 are automatically reproduced by diagonalizing 
the local /-electron term Hi, once we fix the CEF parameters 
in the form of one-electron potential Eq. (7). 
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Fig. 1 . Energies of / electrons as functions of x for (a) the LS coupling 
and (b) the j-j coupling schemes for n=2. The magnitude of the CEF po- 
tential energy is fixed as \W\/U=0.00l, 



2.2 Comparison with LS and j-j coupling schemes 

Before proceeding to the discussion on the energy levels of 
Th point group for n=l~13, we compare the electronic states 
of Hi with those of LS and j-j coupling schemes, since it is 
quite instructive to understand the meanings of the CEF po- 
tential in /-electron systems. We introduce "U" as an energy 
scale for the Racah parameters, A, 5, C, and D. In this sec- 
tion, U is the energy unit, which is typically considered to be 
1 eV. 

In general, in /-electron systems, the magnitude of the 
CEF potential is smaller than both spin-orbit coupling and 
Coulomb interactions. Thus, it is reasonable to consider that 
W is always much smaller than A and U. However, there oc- 
cur two situations, depending on the order to take the limits 
of X/W—hxi and U/W^oo. When the limit of U/W—hxi is 
first taken and then we consider the effect of the spin-orbital 
coupling A, we arrive at the LS coupling scheme. On the other 
hand, it is also possible to take first the infinite limit of X/W. 
After that, we include the effect of Coulomb interaction, lead- 
ing to the j-j coupling scheme. In the present local /-electron 
term Hi, it is easy to consider two typical situations for /- 
electron problems, | Vl^|<gcA<C/ and |W|<C£7<A, correspond- 
ing to the LS and j-j coupling schemes, respectively. 

Let us consider here the case of n=2 as a typical example 
of the comparison with the two schemes. In the LS coupling 
scheme for the / 2 -electron system, we obtain the ground-state 
level as 3 5 with S=l and L=5 from the Hund's rules, where 
S and L denote sums of each /-electron spin and angular mo- 
mentum, respectively. Upon further including the spin-orbit 
interaction, the ground state is specified by J=4 expressed as 
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3 H4 in the traditional notation. Note that the total angular mo- 
mentum J is given by J=\L — S\ and J=L + S for n<l and 
n>7, respectively. In order to consider further the CEF ef- 
fect, by consulting with the table of Hutchings for the case 
of J=4, we easily obtain the nine eigen values, including Ti 
singlet, doublet, and two kinds of T4 triplets, as shown 
in Fig. 1(a). We note that in the LS coupling scheme, we set 
F(4)=60, F(6)=1260, F*(6)=30, y=0.05, and W/U=Q.00l. 
Note also that two T4 triplets in Tj, are obtained by the mix- 
tures of T± and Y§ triplets in the Oh point group. 

In the j-j coupling scheme, on the other hand, first we take 
the infinite limit of A. Thus, we consider only the j=5/2 sex- 
tet, where j denotes the total angular momentum of one /- 
electron. In the / 2 -electron system, two electrons are accom- 
modated in the sextet, leading to fifteen eigen states including 
J=4 nontet, J=2 quintet, and J=0 singlet. Due to the effect of 
Hund's rule coupling, J=4 nontet becomes the ground state. 
When we further include the CEF potential, it is necessary 
to reconsider the accommodations of two electrons in the f 1 - 
electron potential with r> doublet and r§ quartet. Note that 
the difference between Oh and Ti, does not appear for the case 
of J=5/2, namely in the j-j coupling scheme. Thus, in the j-j 
coupling schemes, except for the energy scale W, only rele- 
vant CEF parameter is x, leading to the level splitting between 
Tj doublet and r§ quartet. For the j-j coupling scheme, we 
set F(A)=60 and W/U=-0.001. Note that the minus sign in 
W is added for the purpose of easy comparison with the LS 
coupling scheme. As shown in Fig. 1(b), the J=4 nontet is 
split into Ti singlet, doublet, T4 triplet, and r$ triplet. For 
simplicity, we use the notations for irreducible representations 
of Oh symmetry. The ground state for x>0 is r$ triplet com- 
posed of a couple of Tg electrons, while for a;<0, it is Ti 
singlet which is mainly composed of two ^^ electrons. Note 
that for x>0, the first excited state is triplet, composed of 
Ti and T 8 electrons. 

One may complain that the energy levels in the j-j coupling 
scheme do not agree with those of the LS coupling scheme 
at the first glance. In order to reply to this complaint, let us 
directly diagonalize Hf by changing U and A. Here it is con- 
venient to introduce a new parameter to connect the LS and 
j-j coupling schemes as 

X/\W\ 



'" u/\w\ + x/\w\' (10) 

where both U and A are very large compared with \ W\, since 
we are considering the actual /-electron compound. Note that 
fc=0 and 1 are corresponding to the limits of X/U=0 and 
X/U=oo, respectively. Then, we can control the change of two 
schemes by one parameter k. 

In Figs. 2(a)-(d), we show the energy levels of Hf for sev- 
eral values of k with both A and U larger than \W\. Racah 
parameters are set as A/U=10, B/U=0.3, C/U=0.1, and 
D/U=0.05 in the unit of U. As described above, the CEF 
potential is always small and here we set W/U=— 0.001. In 
Fig. 2(a), results for fc=0.1 are shown. In this case, X/U=Q.l 1 
and the condition A/|W|^1 is still satisfied. Without the 
spin-orbit interaction, the ground-state level is expressed as 
3 H with 5=1 and L=5 due to Hund's rules. When we increase 
A, the multiplet labeled by J is well separated and the ground- 
state level is specified by J=4, as expected from the LS cou- 
pling scheme. Then, the energy levels in Fig. 2(a) are quite 



0.04 



o 

£ 0.02 

§> 0.00 
c 

CD 

g -0.02 

a> 

Lu 



-0.04 



z> 

M — 

O 



0.02 - 



c 

Z5 

>5 0.00 

k_ 

(D 

c 

<D 

c 

(D 

a> 
Lu 



-0.02 , 



-0.04 



° 0.02 

H — 1 

"c 

^ 0.00 

(D 

-0.02 
Lu 



■0.04 



1 1 r- 

(a) k=u.i 


^\ r 












(b) k=0.8 


^ 








r (2) , 


- — r; 






(c) k=0.99 

\ r (i) 


1 4 










r 






^23__ 

\ 1 4 




1 4 



-1.0 



1.00 



/\ 0.95 



\/ 0.90 



0.85 



LS coupling scheme 




1.0 



|W|/U=-0.001 
x=-0.1 
y=0.3 



j-j coupling scheme 



0.0 0.2 0.4 0.6 0.8 1.0 

k 

Fig. 2. Eigen energies of H{ as functions of x for (a) fe=0.1, (b) fc=0.8, and 
(c) fc=0.99. Racah parameters are set as A/U=\Q, 5/(7=0.3, C/U=0.1, 
and D/U=0.Q5. CEF potentials are given by W/U=-0.001 and y=0.3. 
Solid squares at fe=0 and 1 are obtained separately from the LS and j-j 
coupling schemes, respectively. 



similar to those of Fig. 1(a), since we are now in the region 
where the LS coupling scheme is appropriate. 

Even when A is further increased and k is equal to 0.5, the 
structure of the energy levels is almost the same as that of the 
LS coupling scheme (not shown here). However, when k be- 
comes 0.8, as shown in Fig. 2(b), the energy level structure 
is found to be deviated from that of the LS coupling scheme. 
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Rather, it becomes similar to the energy level structure of the 
j-j coupling scheme. To see the agreement with the j-j cou- 
pling scheme, we consider very large A to give k=0.99. As 
shown in Fig. 2(d), we can observe the energy level structure 
similar to Fig. 1(b). Especially, the region of the T 3 ground 
state becomes very narrow, as discussed later. Thus, it is con- 
cluded that Hf correctly reproduces the energy levels both for 
the LS and j-j coupling schemes. We also stress that Hf pro- 
vides correct results for any value of /-electron number. 

A crucial point is that the structure of energy levels is con- 
tinuously changed, as long as A and U are large compared 
with the CEF potential. Namely, the states both in the LS and 
j-j coupling schemes are continuously connected in the pa- 
rameter space. Thus, depending on the situation to consider 
the problem, we are allowed to use the L S or j-j coupling 
scheme. In order to clarify this point, we evaluate the overlap 
(*|*ls). where and |*ls) are me eigenstate of Hf and 
that in the LS coupling scheme, respectively. In Fig. 2(d), we 
show the overlap for the case of T x ground state for x=— 0.1, 
y=0.3, and W/U=-0.001. For fc=0, (*|* LS )=1 due to the 
definition. The overlap is gradually decreased with the in- 
crease of k, but it smoothly converges to the value at k=l, 
i.e., the j-j coupling scheme. Note that the overlap between 
the eigenstates of the LS and j-j coupling schemes is as large 
as 0.865, which seems to be larger than readers may naively 
anticipated from the clear difference between Figs. 1(a) and 
(b). It is not surprising, if we are based on the principle of 
adiabatic continuation, since the eigenstates of the LS and j- 
j coupling schemes are continuously connected. 

Remark that we can observe the common structure around 
at the value of x, in which singlet and triplet ground states 
are interchanged. Namely, essential point of the singlet-triplet 
crossing can be captured both in the two schemes, which will 
be important in the following discussion. However, the r 3 
non-Kramers doublet cannot be the ground state in the j-j 
coupling scheme, since the doublet in the J=4 nontet is com- 
posed of degenerate two singlets formed by T*; and Tg elec- 
trons. As easily understood, such singlets are energetically 
penalized by the Hund's rule interaction and the energy for 
triplet composed of Tr and T$ electrons is always lower 
than that of the singlets. Thus, in the j-j coupling scheme, r 3 
non-Kramers doublet does not appear as the ground state ex- 
cept for x=0. Of course, if j=7/2 octet is explicitly included 
and A is kept finite, it is possible to reproduce T 3 doublet. 
Namely, taking account of the effect of j=7/2 octet is equiva- 
lent to consider the local /-electron term Hf, as we have done 
in this subsection. However, if we expand the Hilbelt space so 
as to include both j=5/2 sextet and j=H2 octet, we lose the 
advantage of the j-j coupling scheme considering only j=5/2 
sextet. 

One may claim that it is possible to reproduce the result of 
the LS coupling scheme even within the j-j coupling scheme, 
just by considering that the CEF potential for J=4 in the LS 
coupling scheme also works on the J=4 / 2 -states composed 
of a couple of / electrons among ,7=5/2 sextet. However, such 
a procedure is not allowed. The reasons are as follows. First it 
should be noted that the CEF potential is not determined only 
by the value of J. For instance, as we will see later, the results 
of the energy levels for n=l and 13 are apparently different, 
even though both of the ground-state multiplets are character- 
ized by J=7/2, since the CEF potential depends also on the 
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values of L and S. Note that for n=7, S=7/2 and L=0, while 
for n=13, L=3 and 5=1/2. For the case of n=2, even if the 
/ 2 -state is characterized by J=4 in the j-j coupling scheme, 
we cannot simply validate the application of the CEF poten- 
tial in the LS coupling scheme to the J=4 / 2 -state in the j-j 
coupling scheme. 

Second we should note again that the CEF effect appears 
only as a one-electron potential. The CEF potential working 
on the two-electron state should be given by the superimposi- 
tion of the one-electron potential. Thus, when we use the basis 
which diagonalizes the spin-orbit interaction, it is necessary 
to consider that the CEF potential should work on the state la- 
beled by the z-component of j. This is the only way to define 
the CEF potential in the j-j coupling scheme, even though the 
T 3 non-Kramers doublet is not reproduced. As mentioned in 
the above paragraph, if j=7/2 octet is included in addition to 
j=5/2 sextet, it is possible to reproduce the results of the LS 
coupling scheme including the non-Kramers doublet. 

2.3 CEF schemes for filled skutterudites 

In the previous subsection, we have examined the local /- 
electronic states in comparison with the results of the LS and 
j-j coupling schemes for n=2. In this subsection, we change 
the /-electron number from n=l to 13, in order to obtain the 
local /-electron states corresponding from Ce- to Yb-based 
filled skutterudites. Note that the results in this subsection 
will be comparable with the experimental ones in the high- 
temperature region, which is equivalent to the atomic limit. In 
the calculation, we set the parameters as A/U=\Q, B/U=0.3, 
C/U=0.l, D/U=0.05, X/U=0A, and W/U=-0.001 in the 
energy unit of U. As already mentioned, we may consider 
C/^leV, indicating that A~0.4eV and | VF|~4meV, which are 
realistic values for rare-earth ions. One of CEF parameters y 
is also fixed as y=Q3 and we evaluate the eigen energies of 
Hf as functions of x. Note also that for each n, the energy is 
appropriately shifted for the easy comparison. 

In Fig. 3, we show the results for n=\, 2, and 3, correspond- 
ing to Ce 3+ , Pr 3+ , and Nd 3+ ions, respectively. For n=l, as 
shown in Fig. 3(a), we observe two curves, T 5 doublet and r 6 7 
quartet. Note that T$ and Tqi in Th correspond to T7 and Tg 
in Oh, respectively. 8 In order to determine x for filled skut- 
terudites, let us move to the case of n=2. In Fig. 3(b), for 
n=2, we can observe four eigen states, as described above. 
Since it has been confirmed experimentally that the ground 

(2) 

state is T\ and the first excited state is T\ with small exci- 
tation energy, it is quite reasonable to set x around at 0.37. 
Within a simple point-charge model, it is possible to fix the 
CEF potential irrespective of n. Thus, when we move back 
to the case of n=\, the ground state is quartet. When we 
turn our attention to Fig. 3(c) for the result of n=3, in the re- 
gion of x>0, the ground state is r 6 7 quartet, consistent with 
the experimental results for elastic constant 36 and magnetic 
entropy 37 of NdFe4Pi2. 

Now let us see the results for n=4, 5, and 6. The case of n=4 
corresponds to Pm 3+ , but Pm is an artificially prepared ele- 
ment. Thus, we show the result only from the theoretical in- 
terest. Since the multiplet including the ground state is charac- 
terized by J=4 for the case of n=4, as shown in Fig. 4(a), the 
structure of the eigen energies is essentially the same as that 
for n=2. Note that the results of x<0 for n=2 appears in the 
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region of x>0 for n=4. Then, the ground state for n=4 should 
be Ti singlet for x>0. The case of n=5 is corresponding to 
Sm 3+ . As shown in Fig. 4(b), the ground state is 1^7 quar- 
tet for x>0, consistent with the experimental results of elastic 
constant 38 and magnetic entropy 39 for S111RU4P12. For n=6, 
the total angular momentum J becomes zero, since S=L=3 
in this case. Thus, the singlet ground state appears without 
the effect of CEF potential. The excitation energy is in the 
order of 0.1, which will be about 1000K in the present en- 
ergy unit. The case of n=6 corresponds to trivalent Eu 3+ , but 
in the filled skutterudite structure, significant contributions 
from Eu 2+ have been found in the measurement of magnetic 
susceptibility 40 and 151 Eu Mossbauer experiment. 41 Divalent 
Eu 2+ ion has / 7 configuration, which is discussed in the fol- 
lowing. 

At half-filling (n=7) corresponding to Eu 2+ or Gd 3+ , the 
result becomes quite simple. Due to the Hund's rule coupling, 
the total spin S is equal to 7/2, while the angular momen- 
tum L is zero. Then, the total angular momentum J is given 
by J=S=l/2, but no effect of CEF potential occurs, leading 
to the J=7/2 octet which is fiat irrespective of x, as shown 
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Fig. 4. Eigen energies of Hf as functions of x for (a) n=4, (b) n=5, (c) 
n=6, and (d) n=l. Other parameters are the same as those in Fig. 3. 



in Fig. 4(d). As is well known in the theory for CEF, since 
the CEF potential works only on the electron charge distribu- 
tion, the case of L=0 is not affected from the outset. Note that 
the excitation energy in the present calculation is estimated as 
5eV. 

Now we will move to the situation for n>l, correspond- 
ing to heavy lanthanide ions. In Fig. 5, we show the results 
for n=8, 9, and 10. The case of n=8 corresponds to Tb 3+ . As 
shown in Fig. 5(a), for n=8, the ground state is Ti singlet 
around at x~0.37. However, experimental results on mag- 
netic susceptibility for TbRu4Pi2 have suggested two succes- 
sive magnetic transitions at 20K and 10K. 20 As long as we 
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=9 and (c) ^' Eigen energies of Hf as functions of x for (a) n=l 1, (b) n=12, and 

(c) n=13. Other parameters are the same as those in Fig. 3. 



fix 2=0.37, the present CEF calculation indicates the Ti sin- 
glet ground state. To explain the discrepancy, for instance, we 
may consider the deviation from the ideal point-charge pic- 
ture. In particular, due to the difference in the ion radius, the 
CEF potential can be different from material to material and 
the values of CEF parameters may be changed. If we slightly 
increase the value of x, the ground state turns to be T4 triplet, 
leading to the magnetic properties at low temperatures ob- 
served in Tb-based filled skutterudites. 

In Fig. 5(b), the result for n=9 corresponding to Dy 3+ ion 
is shown. Around at x=0.37, the ground state is T$ doublet, 
which is T-j doublet in the term of Oh symmetry. Thus, the 
ground state should be magnetic, but the effect of quadrupole 
may appear in Dy-based filled skutterudites, if the CEF po- 
tential is deviated from the present one. In Fig. 5(c), we show 
the result of n=10, corresponding to Ho 3+ ion. In this situa- 
tion, around at a;~0.37, the ground state is triplet, while 
non-Kramers doublet is the first excited state with very 
small energy difference, less than 10K in the present energy 
unit. In addition, as mentioned above, x may be changed from 
0.37, which is the value determined for Pr-based filled skut- 
terudite. Thus, the ground state of Ho-based filled skutteru- 



dite is believed to have magnetic nature, but the effect of 
quadrupole degree of freedom of the doublet may appear 
as an anomaly in physical quantities such as elastic constant. 

Finally in Fig. 6, we show the results for n=l 1, 12, and 13, 
corresponding to Er 3+ , Tm 3+ , and Yb 3+ ions. As observed 
in Figs. 6(a)-(c), the ground state is T$ Kramers doublet for 
n=ll or 13, while Ti singlet for n=12 around at x~0.37. 
Note that these states are not drastically changed in the re- 
gion of 0<a;<0.8, indicating the stable ground states for the 
filled skutterudites with n=l 1 , 12, and 13. 

3. Multiorbital Anderson Model 

In the previous section, we have discussed in detail the local 
/-electron states for n=l~13. However, in order to consider 
the low-energy state, it is necessary to include the itinerant 
effect of / electron. Although the actual material should be 
described by the periodic system, it is a hard task to treat such 
a periodic model by keeping seven /-orbitals. Thus, in this 
section, we consider the impurity Anderson model, by taking 
into account of the hybridization between / and conduction 
electron states. The symmetry of the conduction band is im- 
portant to understand correctly the low-energy properties of 
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filled skutterudite compounds. 

3.1 Hamiltonian 

The Anderson Hamiltonian is written by 



H = 53 £ *> C L C ^+ ^{Vmcljma +h.C.) + H{, (11) 



ho 



where is the dispersion of conduction electron, Cfc CT is the 
annihilation operator for conduction electron with momentum 
k and spin a, and V m is the hybridization between conduction 
and / electrons. The local /-electron term Hf is already given 
in the previous section. In the filled skutterudites, the conduc- 
tion band is given by a u , constructed from p-orbitals of pnic- 
togen. 28 Note that the hybridization occurs between the states 
with the same symmetry. Since the a u conduction band has 
xyz symmetry, we set 



V for m = 2, 
V m = I -V for m = -2, 
otherwise. 



(12) 



Throughout this section, the energy unit is taken as Do, which 
is half of the bandwidth of the conduction band. From the 
band-structure calculation, the bandwidth is estimated as 2.7 
eV in PrRu 4 Pi2, 17 indicating that D =l-35 eV. Note that in 
order to adjust the local /-electron number, we change the 
chemical potential for each n. 

3.2 Method 

In this paper, the Anderson model is analyzed by using the 
numerical renormalization group (NRG) technique. 42 43 In the 
NRG calculations, in order to consider efficiently the conduc- 
tion electrons near the Fermi energy, the momentum space is 
logarithmically discretized and the conduction electron states 
are characterized by "shell" labeled by N. The shell of N=Q 
denotes an impurity site including / electrons. The Hamilto- 
nian is transformed into the recursion form as 

H N+1 = VAH N + t N }^(d N(r c N+1(7 + c' N+1(T c N<T ), (13) 

a 

where A is a parameter for logarithmic discretization, cno- 
denotes the annihilation operator of conduction electron in 
the iV-shell, and t jv is "hopping" of electron between N- and 
(N + l)-shells, given by 

(1 + A- 1 )(1- A-^- 1 ) 



(14) 



2^(1 - A- 2W ^ 1 )(1 - A" 2JV -3) 
The initial term Hq is given by 

H = AT 1/2 [Hi + VM>*fmv + fl a C 0a )]. (15) 
ma 

In this paper, A is set as 5 and we keep 4000 low-energy states 
for each renormalization step. 

In order to see magnetic properties, we evaluate the mag- 
netic susceptibility of /-electron, defined by 

~TrM? AT e~ H "/ T TiS 2 z N er B **/ T ~ 



y= — lim 



Tre- HN / T 



Tj:e -H° N /T 



(16) 



l)/2 



where T is a logarithmic temperature given by T=A ( N 
in the NRG calculation, is the Hamiltonian without Hf, 
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Fig. 7. (a) Excitation energy vs. x around at ie~0.37 for n=2. (b) Magnetic 
susceptibility and (c) entropy of / electron vs. T for n=2 and £=0.36, 
0.365, 0.366, and 0.37. 



and magnetic moment is defined by 

M X>N = -fiB ^2(m + g s a/2)f} na f ma + S Z , N . 

m,(T 

Here is the Bohr magneton, g s =2, and 

JV 

Sz,N = -9sm ^i^/^^naCna- 
n—0 a 

The free energy F for / electron is evaluated by 



F = -T lim 

N~*oo 



In Tre - ^ /T - In Tre~ H « /T 
The entropy S is obtained by S=-dF/dT. 



(17) 



(18) 



(19) 
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3.3 Result for n=2 

First let us briefly review the results for n=2, 44 ' 45 corre- 
sponding to the Pr-based filled skutterudites. As mentioned 
in the previous section, it has been confirmed experimentally 

(2) 

that the ground state is Fi and the first excited state is T\ 
with small excitation energy for Pr-based filled skutterudites. 
In Fig. 7(a), we show the excitation energy for n=2 around at 
2=0.37 in a magnified scale. Then, we consider the region of 
0.36<2<0.37 in the following. 

In Fig. 7(b), we show calculated results of T\ for n=2 and 
several values of x between 0.36 and 0.37. For x<0.365, T\ 
becomes zero for small T, while for x>0.366, T\ takes con- 
stant at low temperatures. 44 As shown in Fig. 7(c) for the 
calculated results of entropy S, we obtain log 3 as residual 
entropy for 2=0.366 and 0.37, indicating that local triplet re- 
mains in the low-temperature region. Of course, the local mo- 
ment is eventually suppressed at low temperatures in actuality. 
For instance, there should exist the hybridization with con- 
duction bands other than a u . What we emphasize here is that 
magnetic fluctuations still remain at relatively low tempera- 

(2) 

tures even when T± is the ground state, if T\ is the excited 
state with small excitation energy. In the following calcula- 
tion, we set 2=0.366. 

We remark that the transition between magnetic and non- 
magnetic states is governed by the exchange interaction J c f 
between / electrons and the conduction band, expressed as 
J c f = V 2 /SE, where SE denotes the energy difference be- 
tween the f 2 and / 3 (or f 1 ) lowest-energy states. In fact, it 
has been found that the boundary curve between magnetic and 
non-magnetic phases is proportional to V 2 , 45 suggesting that 
the dominant energy scale should be J c f for the appearance of 
magnetic fluctuations. Note that the transition between mag- 
netic and non-magnetic states seems to be abrupt, not gradual 
as observed in the usual Kondo system. In the Kondo prob- 
lem, the local moment is suppressed only by hybridization 
with conduction electrons, but in the present case, a singlet 
ground state is also obtained through local level crossing due 
to the CEF potential, in addition to the hybridization process. 
Furthermore, localized orbitals exist in the present model, as 
will be discussed in the next subsection. Namely, the mag- 
netic moments of the / orbitals hybridized with the conduc- 
tion band are suppressed as in the Kondo effect, while the 
moments of the localized orbitals are not screened. An abrupt 
change in the magnetic properties is caused by the duality of 
the / orbitals in combination with the level crossing effect 
due to the CEF potential. 

3.4 Reconsideration in a j-j coupling scheme 

In order to understand the physical meaning of the result 
for n=2, it is useful to express the /-electron state with xyz 
symmetry in the j-j coupling scheme as 

|xyz,<7) - y/A/7\7/2,T 7 ,a) - y^J7\5/2,T 7 ,a), (20) 

where \j, T, a) denotes the state in the j-j coupling scheme 
with total angular momentum j, irreducible representation T 
for the Oh point group, and pseudospin a. We note that in the 
j-j coupling scheme, we use the irreducible representation of 
Oh, since the CEF potential for one /-electron state in the 
j=5/2 sextet is the same for Oh and T^. Note also that there is 
one-to-one correspondence between a and u. The above rela- 



tion indicates that only the T 7 state is hybridized with a u con- 
duction band states with xyz symmetry, leading to the Kondo 
effect, while the r 8 electrons are localized and become the 
source of local fluctuations at low temperatures. It is an im- 
portant issue of filled skutterudite structures that the nature of 
the / electrons can clearly be distinguished as itinerant T 7 or 
localized T 8 . This point provides a possible explanation for 
the heavy-fermion phenomenon occurring in f 2 -electron sys- 
tems such as Pr-based filled skutterudites. 

In the j-j coupling scheme, we generally assume the large 
spin-orbit coupling, indicating that j=7/2 octet is discarded 
and only j=5/2 sextet is taken into account. We also note that 
the conduction band has T 7 symmetry in the filled skutteru- 
dite compounds. Then, the Anderson Hamiltonian in the j-j 
coupling scheme is written as 

kcr ka 

where E k is the dispersion of T 7 conduction electron, the en- 
ergy unit in this subsection is also taken as D which is half of 
the bandwidth of the conduction band, a ka is the annihilation 
operator for T 7 conduction electron with momentum k and 
pseudospin a, /j 7<7 is the annihilation operator for j=5/2 f- 
electron on the impurity site i with spin a and "orbital" 7, V 
is the hybridization between conduction and /-electrons with 
T 7 symmetry, and H\ denotes the local /-electron term at site 
i. Note that the orbital index 7 is introduced to distinguish 
three kinds of the Kramers doublets, two r 8 and one T 7 . Here 
"a" and "b" denote two r 8 and "c" indicates T 7 . 

The local /-electron term Hi in the j-j coupling scheme is 
given as 21 

7,(7 7l~74Cl,0'2 

X /i7i (Ti /i72 <j 2 />73 CT 2 /174 ci 1 (22) 

where B 7 is the CEF potential for j=5/2. Since the CEF po- 
tential is already diagonalized, it is convenient to introduce a 
level splitting between T 7 and T 8 as 

A = Br 8 - B Fr = QWx. (23) 

The Coulomb integral -f 7 r 1 1 j ' 7 r 2 2 i73j74 m the j-j coupling scheme 
is expressed by using other Racah parameters, E , E\, and 

fi 1 21 

First let us discuss the local /-electronic properties of H\, 
when we change the level splitting A, although the essen- 
tial point has been already discussed in Sec. 2. As shown 
in Fig. 8(a), the ground state is T 5 triplet for A<0, while it 
is Ti singlet for A>0. Note that for positive A, the triplet 
excited state is T4, composed of F7 and T 8 electrons. Note 
again that we use the irreducible representation of Oh, since 
there is no difference between Oh and Th in the j-j coupling 
scheme. We again apply the NRG method to the Anderson 
model in the j-j coupling scheme. In Figs. 8(b) and (c), we 
show the results of T\ and the entropy S of / electrons for 
several values of A. Racah parameters are set as Eq/Do=5, 
Ei/Dq=2, and E2/ Dq=0.5, which are appropriate for rare- 
earth compounds. The hybridization is set as V/ Dq=0.05. 
For high temperatures, since all the J=4 states contribute to 
the /-electronic properties, we obtain the entropy of log 9. 
For A/D >6 x 10~ 5 , both susceptibility and entropy even- 
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Fig. 8. (a) Eigen energies of H\ for n=2. Insets denote /-electron config- 
urations for Ti, T4, and Ts states, in which two / electrons are accom- 
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electron vs. T for n=2 based on the Anderson model in the j-j coupling 
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tually go to zero in the low-temperature region, while for 
A/Dq<5 x 10 -5 , magnetic fluctuations remain significant 
at low temperatures. Note that the residual entropy is log 3 for 
A/D <5 x 1CT 5 . 

Now we consider the reason why the essential magnetic 
features of Pr-based filled skutterudites can be captured as 
described in the j-j coupling scheme. As shown in Fig. 8(a), 
when the ground state is a Ti singlet, there are two triplet ex- 
cited states, T4 and IV The Ts triplet is composed of a pair of 
Tg electrons, while the triplet is composed of one and 
one Tg electron. Since Tg electrons do not hybridize with 1^ 
conduction electrons, the T$ triplet can survive. Note that 
triplets in Ti, are given by the mixtures of and T$ in Oh- 
Such a mixing is not included in the j-j coupling scheme, 
but there exists a T$ excited state. Thus, the local triplet still 
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Fig. 9. (a) Magnetic susceptibility and (b) entropy of / electron vs. T for 
n=l~6. Racah parameters are set as A/D =10, B/D =0.3, C/D =0.\, 
and D/Do=0.05. Spin-orbit coupling and of hybridization are set as 
A/Do=0.4 and V/ Dq=0.Q5, respectively. As for CEF parameters, we set 
W/A)=-0.001, z=0.366, and y=0.3. 



remains, as long as the excitation energy is smaller than J c f . 

Summarizing this subsection, the low-temperature /- 
electron properties can be well captured even in the Ander- 
son model constructed from the j-j coupling scheme, in the 
parameter region in which singlet and triplet ground states 
are interchanged. In particular, even though the ground state 
is Ti singlet, the magnetic fluctuations survive in the low- 
temperature region, when the triplet excited state exists with 
very small excitation energy. Thus, the j-j coupling scheme 
is applicable to investigate the microscopic /-electron prop- 
erties of rare-earth filled skutterudites. 

3.5 Results for n=l ~i5 

Thus far, we have focused on the case of n=2, but we can 
further study the cases of n=l~13 based on the original An- 
derson model Eq. (11). In the point-charge picture, the CEF 
parameters do not depend on n, since the CEF effect is due to 
electrostatic potential from ligand ions surrounding the rare 
earth. In the following, we fix x=0.366, in which the local 
f 2 ground state is a Ti singlet, but for which we have found 
significant magnetic fluctuations at low temperatures. 

In Fig. 9(a), magnetic susceptibilities for n<7 are shown 
for 2=0.366. For n=2, 4, and 6, the ground state is Ti singlet, 
indicating that both magnetic and orbital fluctuations should 
be rapidly suppressed with decreasing temperature. However, 
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(2) 

for n=2, since the T\ triplet is existing with very small ex- 
citation energy, magnetic fluctuations significantly remain, as 
mentioned in the previous subsection. Note that for n=A and 
6, the excitation energy is large and T\ becomes zero imme- 
diately in the low-temperature region. On the other hand, for 
n=l, 3, and 5, the ground state is Tg7 quartet, as confirmed 
from the residual entropy of log 4, as shown in Fig. 9(b). Here 
we remark that the case of n=5 corresponds to Sm 3+ ion. In 
NMR experiments for Sm-based filled skutterudites, 46 it has 
been suggested that magnetic fluctuations seem to appear in 
the low-temperature region, consistent with the present nu- 
merical result. 

Next we turn our attention to the cases of half-filling (n=7) 
and more than half-filling (n>7). In Fig. 10(a), we show the 
NRG results for susceptibility. First of all, the absolute values 
of x are much larger than those for n<7, since the magnetic 
moment becomes large for n>l due to the Hund's rule in- 
teractions. Typically, at half-filling, total spin S(=J) is equal 
to 7/2 and the Curie constant for an isolated ion is as large 
as 21 /ig/fce- In the wide temperature region, this value has 
been observed in Fig. 10(a) for n=7, indicating that 5=7/2 
spin survives at relatively low temperatures. In fact, we clearly 
see the entropy of log 8 from the octet of 5=7/2, as shown in 
Fig. 10(b). In the extremely low- temperature region, however, 
eventually 5=7/2 spin is screened and only 5=1/2 remains, 
leading to the residual entropy log 2. 



For the cases of n=8 and 12, as discussed in the previous 
subsection, the ground state is T% singlet and thus, the sus- 
ceptibility rapidly goes to zero. Note that in the case of n=8, 
when we change the value of x, the effect of the magnetic ex- 
cited state becomes significant, as in the case of n=2. In order 
to focus on the case of Tb-based filled skutterudites, it is nec- 
essary to analyze in detail the experimental results, but it will 
be discussed elsewhere in future. For n=9, 11, and 13, the lo- 
cal ground state is Kramers doublet, which is the mixture 
of Tq and of Oh- Since Tq state does not hybridize with 
the a u conduction band, the magnetic moment from Tq still 
persists even in the low-temperature region. In fact, we ob- 
serve the residual entropy of log 2 in these cases, as shown in 
Fig. 10(b). 

For n=10, the local ground state is T4 triplet, but it just re- 
mains as it is in the low-temperature region, as observed in 
Figs. 10(a) and (b). This is understood as follows: It is con- 
venient to reconsider the states in the j-j coupling scheme 
and Oh symmetry. The Tr, triplet in the Oh symmetry is com- 
posed of a couple of r§ electrons, while the T4 triplet in the 
Oh symmetry is composed of and r§ electrons. Namely, 
the F5 triplet can survive, since Tg electron does not hybridize 
with T7 conduction electron. In the Th symmetry, T4 triplet 
states are always given by the mixture of and triplets 
of Oh and thus, the r$ component still remains, even after the 
hybridization with conduction electrons. 

4. Superconductivity of Pr-Based Filled Skutterudites 

In order to discuss the appearance of superconductivity 
from a microscopic viewpoint, here we simply consider the 
Hubbard-like model for /-electron systems, 21 in which T7 
electrons move around the system through /-/ hopping in 
the present case. Note that it is necessary to treat the peri- 
odic Anderson model in order to consider simultaneously the 
formation of heavy quasi-particle and the occurrence of su- 
perconductivity. However, for the present purpose to discuss 
the nature of superconducting pairing, we believe that the con- 
struction of the Hubbard-like model for /-electron systems is 
useful. The Hamiltonian is given by 



H = -t 



E 



(/icer/jc 



(24) 



where t denotes an effective hopping amplitude of elec- 
tron between nearest neighbor sites, (i, j), in the bcc lattice. 
In this section, t is taken as a new energy unit. We believe 
that this Hamiltonian can be a canonical model to investigate 
superconductivity of filled skutterudites from the microscopic 
viewpoint. 

Now we evaluate superconducting pair susceptibility P. 
For the purpose, we define the pair operator $ as 



$ = 



i7<ri'7 / 'a' 



^i-yai' ^7 ' <j' ^i7fri / 7 / cr / 1 



(25) 



where $i 7cri / 7 / (T / = }i llT h'Y^ and Vi7<ri'7'<r' is the coeffi- 
cient determined by the eigen state of the maximum eigen 
value P max of the matrix P, given by 

f 1/T - M 
Pi^i'V^'JAadW = / dr < $ i7<Ti'7'<T'( T )*Ui'uV>- ( 26 ) 



Here A(r)=e HT Ae Ht for an operator A. Information on the 
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Fig. 11. (a) Eigen values of pair susceptibility vs. A in the 2-site model for 
T/i=10 -3 . (b) Eigen values of P as functions of magnetic field A/t=0.01 
in the 2-site model, (c) Phase diagram for superconducting pair in the 2- 
site model. Solid curve indicate the boundary between singlet and triplet 
pairs, while the broken curve denotes the boundary between two singlet 
pair states. 



symmetry of Cooper pair is included in (p. 

In Fig. 11(a), eigen values of pair susceptibility are plot- 
ted as functions of A for the 2-site case. Since three or- 
bitals are included at each site, it is difficult to enlarge the 
system size for the evaluation of susceptibility, but we can 
grasp some essential points of pair wavefunction in an un- 
biased manner. First we observe that only in the region of 
A/i«0, P max is significantly enhanced. The corresponding 
eigen state indicates an off-site singlet pair of Tg electrons. 
For A/i<§Cl and A/0>1, we obtain off-site and on-site sin- 
glet pair of electrons, respectively. Here we emphasize 
a concept that the appearance of superconductivity is con- 
trolled by A. When | A| is increased, pair susceptibility is sup- 
pressed in any case, since insulating states are dominant. For 
positive A with A/t^>l, T7 is doubly occupied, leading to 
charge-density-wave (CDW) state, while for negative A with 
|A|/£^>1, spin 5=1 is formed by a couple of r§ electrons, 
leading to spin-density-wave (SDW) state. When two differ- 
ent insulating phases compete with each other, a metallic state 
appears in the competing region 47 and the superconductivity 
is expected to appear there. 45 

Let us see what happens under the magnetic field. We eval- 
uate the pair susceptibility by adding the Zeeman term to the 



Hubbard model, given by 



(27) 



!■/< 



where gj is the Lande's g-factor (gj=6/7 for j=5/2), H 
is an applied magnetic field, /i is the z-component of to- 
tal angular momentum j=5/2, and is the annihilation 
operator of /-electron labeled by /1 at site i. For conve- 
nience, we introduce the non-dimensional magnetic field as 
h=gjfiBH/ (kst). Note that in actual calculations, we use the 
/-electron basis so as to diagonalized the CEF term. As shown 
in Fig. 11(b), the state of -Pmax at A/<=0.01 is found to have 
odd-parity triplet pair in the region of a small magnetic field. 
When we increase the magnetic field, there occurs another 
singlet off-site pair state composed of T? and Tg electrons. 

Based on the 2-site calculation results, we depict the phase 
diagram for superconducting pair state in the plane of A and 
magnetic field H, as shown in Fig. 1 1(c). Since we are inter- 
ested only in the possible symmetry of Cooper pair, we do not 
consider the competition among superconducting state and 
other ordered phases. For instance, when the magnetic field 
becomes large enough, quadrupole ordering may be domi- 
nated, as deduced from the experimental result. In the re- 
gion of A<0, for zero magnetic field, there appears the sin- 
glet off-site pair state composed of Tj electrons, as shown in 
Fig. 11(a). However, for A/t<— 0.8, the triplet pair state ap- 
pears only by applying a small magnetic field. Since this state 
is considered to be dominated by the magnetic phase, super- 
conductivity may disappear in actuality. In fact, the maximum 
eigen value in this region is not enhanced in comparison with 
that in the region of A«0. For A>0, we observe a wide region 
of singlet state, which will be dominated the CDW-like phase. 
The local Ti singlet is robust and it is not easily destroyed by 
the magnetic field. Again the maximum eigen value in this re- 
gion is not enhanced and superconductivity is not expected to 
occur, at least, in the present model without phonons. 

Finally, let us focus on the region around A«0. For H=0, 
there occurs the singlet off-site pair, composed of a couple of 
r§ electrons. When we apply a small magnetic field in the de- 
generate region around A«0, we can observe the odd-parity 
triplet pair. The maximum eigen value in this region is clearly 
enhanced in comparison with that for A>0 or A<0 and thus, 
it may be connected to the superconductivity in the thermo- 
dynamic limit. Upon further increasing H, the triplet phase 
eventually disappears and a singlet-pair state turns to occur, as 
mentioned in Fig. 1 1(b). However, such a state will be domi- 
nated by quadrupole ordered phase. 

Readers may consider that the change from singlet off-site 
to odd-parity triplet pair is related to the second transition in 
the superconducting phase of PrOs4Sbi2 detected by thermal 
conductivity. 11 This is an interesting scenario to explain the 
multiple superconducting phase, but the second transition be- 
low T c has not been confirmed by other experimental mea- 
surements. In any case, it seems to be premature to conclude 
the pairing symmetry only from the present 2-site calcula- 
tions. 

5. Discussion and Summary 

In this paper, we have analyzed two kinds of the Anderson 
models with active orbital degrees of freedom by using the 
NRG technique. It has been found that magnetic fluctuations 
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significantly remain, if T 4 triplet is the excited state with 
small excitation energy. By further analyzing the Hubbard- 
like model constructed from the j-j coupling scheme, we 
have proposed a scenario that anisotropic Cooper pair occurs 
in the interchanged region between singlet and triplet states. 

As mentioned repeatedly in this paper, magnetic fluctua- 
tions remain for the degenerate region in which singlet and 
triplet ground states are interchanged. Such magnetic fluc- 
tuations may be antiferromagnetic, since local triplets easily 
lead to antiferromagnetic insulating state in the periodic lat- 
tice. Thus, the off-site singlet pair appearing in the region of 
A«0 in Fig. 1 1(a) may be the d-wave Cooper pair mediated 
by magnetic fluctuations. Note that in the exciton mechanism 
for the superconductivity of PrOs4Sbi2, d-wave pair was also 
suggested. 50 

However, it is premature to conclude the d-wave pair, since 
the present calculation has been done only in the 2-site case. 
In addition, we are considering the pair in the degenerate re- 
gion, where spin and orbital fluctuations should play some 
roles. Note that T 8 electrons have spin and orbital degrees of 
freedom. 21 In particular, it has been pointed out that triplet 
pair can be induced by the cooperation between spin and or- 
bital fluctuations. 48 ' 49 Within a random phase approximation, 
it is easy to show that the spin and orbital fluctuation terms 
compete each other for spin-singlet channel, while for spin 
triplet pairing, orbital fluctuations are cooperative with spin 
fluctuations. Thus, it is considered that the spin-triplet pair 
occurs in the orbital degenerate systems. In fact, by apply- 
ing a small magnetic field, we have observed the odd-parity 
triplet pair, which may be mediated by cooperative fluctua- 
tions of spin and orbital degrees of freedom. In order to see 
whether the triplet pair is stabilized or not for zero magnetic 
field, it is necessary to solve the Eliashberg equation in the 
thermodynamic limit. This is one of future tasks. 

We should note that the quadrupole component is also in- 
eluded in the T\ ' triplet, although we have emphasized mag- 
netic fluctuations in this paper. In fact, an important role 
of quadrupole fluctuations has been suggested experimen- 
tally 7 ' 51 and theoretically. 52 In the future discussion on su- 
perconductivity, it is necessary to clarify the different roles of 
magnetic and quadrupole fluctuations. 

In the present scenario, unconventional superconductivity 
should disappear when A becomes large. For large A(>0), 
charge degree of freedom is dominant and it should be cou- 
pled to phonons. If we further include electron-phonon in- 
teraction into the present model, conventional s-wave super- 
conductivity may occur for large A. Namely, when A is in- 
creased, unconventional superconducting state is eventually 
changed to conventional one. Such a change has been exper- 
imentally observed in Pr(Osi_ x Ru x )4Sbi2- 15 When x is in- 
creased from to 1, the excitation energy between T± and 
Y\ continuously increases and the superconducting proper- 
ties are changed from unconventional to conventional, consis- 
tent with our scenario. It is another future problem to repro- 
duce quantitatively the change of superconductivity based on 
the Hubbard model Eq. (24) with electron-phonon interaction. 
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